GENERALIZED HUNTER-SAXTON EQUATION AND 
THE GEOMETRY OF THE GROUP OF CIRCLE 
DIFFEOMORPHISMS 
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Abstract. We study an equation lying 'mid-way' between the periodic Hunter- 
Saxton and Camassa-Holm equations, and which describes evolution of rotators in 
liquid crystals with external magnetic field and self-interaction. We prove that it 
is an Euler equation on the diffeomorphism group of the circle corresponding to a 
natural right-invariant Sobolev metric. We show that the equation is bihamiltonian 
and admits both cusped as well as smooth traveling-wave solutions which are natu- 
ral candidates for solitons. We also prove that it is locally well-posed and establish 
results on the lifespan of its solutions. Throughout the paper we argue that despite 
similarities to the KdV, CH and HS equations, the new equation manifests several 
distinctive features that set it apart from the other three. 
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1. Introduction 
The equation we study in this paper is 
(//HS) - u txx = -2/j,(u)u x + 2u x u xx + uu xxx 

where u = u(t,x) is a time-dependent function on the unit circle S 1 = R/Z and 
ju(tt) = J sl udx denotes its mean. The closest relatives of this equation are the 
Camassa-Holm equation 

(CH) u t - u txx + 3uu x = 2u x u xx + uu xxx , 

and the Hunter-Saxton equation 

(HS) Uf xx — 2u x u X x ""nil 

both of which have recently attracted a lot of attention among the integrable systems 
and the PDE communities. The Camassa-Holm equation, originally derived in |FF| as 
an abstract equation with a bihamiltonian structure, was introduced independently in 
|CH| as a shallow water approximation^. The Hunter-Saxton equation first appeared 
in HSJ as an asymptotic equation for rotators in liquid crystals. We show that (//HS) 
can be viewed as a natural generalization of the rotator equation once we allow 
interactions of rotators and an external magnetic field. 

In order to compare equation Q^HSp with the other two equations we first recall a 
few facts about the latter. It is known that for spatially periodic functions (1CHI) is the 
geodesic equation on the infinite-dimensional group T> s (S 1 ) of orientation-preserving 
diffeomorphisms of the unit circle S 1 ~ R/Z of (sufficiently high) Sobolev class H s 
and endowed with a right-invariant metric given at the identity by the H 1 inner 
product, see [M] . 

(u,v)h 1 = / (uv + u x v x )dx. 



Equation flCH) is a completely integrable system with a bihamiltonian structure and 
hence it possesses an infinite family of commuting Hamiltonian flows, as well as 
an associated sequence of conservation laws, see e.g. |CM| for the periodic case. 
Furthermore, it admits soliton-like solutions (called peakons) in both the periodic 
and the non-periodic setting, see [CH]. 

Equation (|HS|) . on the other hand, describes the geodesic flow on the homogeneous 
space of the group V s (S 1 ) modulo the subgroup of rigid rotations Rot(S 1 ) ~ S 1 
equipped with the H 1 right-invariant metric, which at the identity is 



{u,v)^ = \ u x v x dx. 



This equation possesses a bihamiltonian structure as well, see e.g. [HZ} IKM] , 

In Section [2] we will show that the //HS equation has a similar geometric origin: 
it describes the geodesic flow on T> S (S 1 ) with the right-invariant metric given at the 
identity by the inner product 

(u, v) = /j,(u)fj,(v) + / u x v x dx. 



(see also Section [5]). One motivation for considering this inner product comes from 
the observation made in [L2] that there is a natural isometry between T> S (S 1 )/ S 1 
equipped with the H 1 metric and an open subset U of the (infinite-dimensional) 
unit sphere in L 2 (S l ). Consequently, the group V s (S 1 ) can be viewed as an S 1 - 
bundle over U with (•,•) inducing the simplest possible metric in this bundle which 



^For an interesting discussion of the attendant physical principles see [J] . 
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projects to the H 1 metric on the base U. (In a sense, this picture can be thought 
of as an infinite-dimensional version of the Hopf S^-fibration.) Thus equation Q/iHS[ ) 
as a geodesic equation for P S (S' 1 ) with this metric is expected to have the 'nicest' 
properties among the integrable equations related to this group, owing to constant 
sectional curvature of the base U and the 'constant' metric along the fibers S 1 . 

In Section [3] we will present two variational principles for the /iHS equation and 
show how it arises as an asymptotic rotator equation in a liquid crystal with a pre- 
ferred direction if we take into account the reciprocal action of dipoles on themselves. 
The bihamiltonian structure of Q^HS| ) together with its attendant infinite hierarchy 
of conservation laws will be described in Sections U] and [H 

In spite of the presence of the nonlocal expression fj,(u) the /iHS equation is a well- 
defined PDE (rather than a mixed integro-differential equationj^. This is because 
the mean of any solution u is constant in time and hence is completely determined 
by the initial condition. In light of this fact the term 2/i(u)u x may appear rather 
harmless and might be expected to have only a minor influence on the phenomenology 
of the equation. However, in Section [5] we shall see that it has a strong effect on 
well-posedness of the associated Cauchy problem in that it is responsible for Q/iHS[ ) 
admitting global (in time) solutions. In the same section a blow-up mechanism for 
/iHS solutions will also be described. In this respect the /iHS equation resembles 
(|CH|) rather than (|HS|) as all classical solutions of the periodic HS equation except 
space-independent solutions blow up in finite time, see [Gl] or [Y]. Let us point out 
that just like the breaking waves of (|CH|) can be continued as global weak solutions 
(see |BC2| ) the solutions of (|HSp can also be extended to all times (see [BC| or |L3| 
for the periodic case). 

Furthermore, whereas (|HSp does not have any bounded traveling-wave solutions 
at all, the /iHS equation will be shown in Section [6] to admit traveling waves (in the 
periodic case) that can be regarded as the appropriate candidates for solitonsO 

In Section [7] we compute the sectional curvature of V s (S 1 ) and in Section [8] we 
revisit the derivation of the /iHS equation on the Virasoro group and describe its two 
compatible Poisson structures. 

Finally, we remark that although our focus in this paper is on the periodic case 
many of the arguments presented here can be adapted, with obvious modifications, 
to study the corresponding equation on the real line. 



2. EULERIAN NATURE OF THE /iHS EQUATION 

2.1. Geodesic flow on the group of circle diffeomorphisms. Let V s (S 1 ) de- 
note the set of circle diffeomorphisms of a sufficiently high Sobolev class H s and let 
TV'iS 1 ) be its tangent bundle^ Consider the following inner product on the cor- 
responding Lie algebra of vector fields, that is, on the tangent space TidD s (S l ) = 

2 We mention here that a related (nonperiodic) equation obtained by replacing fi(u) by a constant 
has appeared previously as a special case of a model derived in MN for the propagation of surface 
waves in an ideal fluid with surface tension, and whose soliton solutions were studied using tools of 
algebraic and complex geometry in [ACHM] . 

^That (|HS|) does not admit any bounded traveling wave solutions is related to the fact that the 
equation itself is defined on the homogeneous space T> a (S 1 )/ S 1 , i.e. 'modulo rotations' and hence 
one considers all solutions modulo traveling. 

4 In the bulk of the paper we can assume that the Sobolev index s > 3/2 is simply large enough 
for the computations to make sense (for example, taking s = oo, which corresponds to the Frechet 
topology on smooth functions, will do). In Section [5] on well-posedness, where analytical details are 
of essence, we will be more precise. 
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{u{x)d/dx : u G H S {S V )}: 

uAv dx . 



(2.1) (u, v) = h(u)/j,(v) + J u x v x dx 



where the inertia operator A = [i — d x acts on a function v by Av = fj,(v) — v xx . 
Using right translations we extend this inner product to a right-invariant metric on 
the group V s (S 1 ) so that 

(2.2) (U, V) v = [ Ur] x dx [ Vt] x dx + f dx 



for any diffeomorphism n G V s (S 1 ) and any U, V G T V V S (S 1 ). 

Theorem 2.1. The [iHS equation is an Euler equation on T> s (S 1 ). More precisely, 
it corresponds to the equation of the geodesic flow on T> S (S 1 ) with respect to the right- 
invariant metric \2. 



Remark 2.2. Observe that A is invertible as an operator from H s onto H s 2 with 
the inverse A~ 1 w = u given explicitly by 



(2.3) 



u ( x ) = ( y - 1 + y^Vh + ( x - \) j J o w ( s ) dsd y 



y r ry rs 

w(s) dsdy + / w(r) drdsdy. 



io Jo J Jo Jo 

Moreover, note that A -1 commutes with d x and n(u) = ^l(Au). 

Proof of Theorem l2.1i Using (|2.ip and the definition of the coadjoint operator 

(ad^(u),w) = (u, [v,w]), 

where [v, w] = v x w — vw x , we readily compute 

(2.4) ad* v {u) = A~ l {2v x Au + vAu x ). 

It follows that the corresponding Euler equation is 

u t = -ad* u (u) = -A~ 1 (2ii(u)u x - 2u x u xx - uu xxx ), 

or, applying A to both sides, equivalently 

H{utj- u txx = -2fi(u)u x + 2u x u xx + uu xxx . 

Integrating both sides of this equation over the circle and using periodicity implies 
that fi(ut) = 0, thus yielding the /uHS equation. □ 

An alternative derivation proceeds by calculating the first variation of the energy 
functional of the right-invariant metric (|2.2|) 



1 



m = 2 J <tf*)> 

defined on the space of smooth paths t — > £(t) between £(0) = id and £(1) in the 
diffeomorphism group. The resulting geodesic equation in V s (S 1 ) right-translated to 
the tangent space at the identity is the ^HS equation. 

Combining Theorem 12.11 with a local well-posedness result of Section [5] we will 
show that the metric (|2.2p has a well-defined Riemannian exponential map on T> S (S 1 ) 
whenever s > 3/2 (Corollary 15. 3p . 
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2.2. fiHS as an evolution equation. Recall that the periodic Hunter-Saxton equa- 
tion 

Htxx ^^x^xx ^VUxxx — 

is defined only as an evolution equation on the quotient T >S (S 1 )/S 1 , i.e. on the space 
of cosets {u(x) + c | x £ S , c £ R}, see e.g. [KMj . The ambiguity introduced by 
arbitrariness of c disappears in the case of the (iHS equation. Indeed, integrating 
equation ( |/iHS[ ) in x gives 

(2.5) u tx = 2fj,(u)u - uu xx - ^u 2 x + r (t), 

where ro(t) is uniquely defined at each time t by the condition that the right-hand 
side be a total x-derivative. Explicitly, one has ro(t) = — 2 / u(u) 2 — ^fi(u 2 ). Note that 
7"o (i) = ro is constant in time, since — 2ro(i) = 4^(u) 2 +/j,(u x ) = 3fi(u) 2 + (u, u), while 
both the //-value and the energy are conserved along ^HS solutions. Integrating once 
again over [0, x] gives 

u t {t,x) = J ^2^(u)u-uu xx -^ul^ dx + r x + q (t). 

The same procedure in the case of the HS equation leads to the arbitrary integration 
constant qo(t), which is the reason for HS being defined only on cosets. This time, the 
condition fi(ut) = implies that the mean value of the expression on the right-hand 
side above must be zero. This defines qo(t) uniquely and hence gives a well-defined 
evolution on the entire space of periodic functions u : S 1 — > IR. 

3. Variational principles and physical motivation 
3.1. fiHS as a bi-variational equation. 

Theorem 3.1. The fiHS equation satisfies two different variational principles. 

We would like to show that (/iHS) arises as the equation 5S = for the action 
functional 



S = J (J Cdx)dt 



with two different densities C. 

Consider the first Lagrangian density 

(3.1) C = u x ut + uu x + 2fi(u)u 2 . 

This is a generalization of the HS Lagrangian density where the /U-term is added. The 
variational principle 5S = gives the Euler-Lagrange equation 

(3.2) ut x = ~\^x ~ uu xx + 2/i(u)u + fi(u 2 ). 

The right-hand-side of this equation has to be understood up to addition of a (possibly 
time-dependent) constant to make it a total derivative: 

(3.3) u tx = --u x - uu xx + 2fi(u)u + /u(n 2 ) + c . 
Here c is defined uniquely by the condition 

J u 2 x dx + 4^(n) 2 + 2fi(u 2 ) + 2c = , 



vanishing of the mean of the right-hand-side in (|3.3p . This immediately leads to 
equation (|2.5[) . Equivalently, by taking the x-derivative, equation (|3.2p can be directly 
understood as the /xHS equation for Uf XX - 
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Remark 3.2. The additional constant term c in equation f|3.3[) is necessitated by the 
term u x ut in the Lagrangian. Indeed, consider any variational principle 



5 (^j j u x utdxdt + f(u) \ = 0. 



It gives the equation ut x = —^df/Su, which makes sense on the circle only if the mean 
of the function 5f jbu is zero. Thus by including the term u x ut in the Lagrangian, 
we are assuming that the Euler-Lagrange equation will describe the evolution of the 
complete derivative u x in the form 

Utx = X(u), 

for some vector field X. By adding an appropriate c for each moment of time to the 
right-hand-side, one can force the result to be a derivative. 

For the Lagrangian (|3.ip this constant will automatically give the appropriate value 
of ro already discussed in Section [2T2J Most importantly, this constant is non-essential 
and will disappear once we differentiate both sides to obtain the equation for ut xx ^ as 
in the original fiHS equation^ 

The second variational representation can be obtained from the Lagrangian density 

(3.4) £ = -{u 2 x + fi(u)u) + ir(z t + uz x ) , 



which is a /x-generalization of the density presented in [HZ] . Varying the correspond- 
ing action with respect to u, it, respectively z, we find 

(3.5) -u xx + n(u) + irz x = 0, 

(3.6) z t + uz x = 0, 

(3.7) -TT t - (ttu) x = 0. 

Thus, using (|3.5h . the expression ut xx + 2u x u xx + uu xxx equals 
ir t z x + irz tx + 2u x (n(u) + irz x ) + u(irz x ) x . 



After employing (13. 6j) and (|3.7j) to replace zt x respectively 7rj, this simplifies to just 
2fj,(u)u x , showing that the equations of motion is the ^HS equation. 

As we will see in Section 14.21 these bivariational principles correspond to the two 
Hamiltonian formulations of the ^HS equation. 



3.2. Physical motivation. First recall the derivation of the HS equation from 
One considers the director field n(x, y, z, t) which describes the average orientation 
of the long rod-like molecules in a nematic liquid crystal. Assuming invariance under 
spatial translations and no fluid flow, the action is given by (see (A. 14) in |HSJ) 

i|n t | 2 - W(n, Vn) + ^\\n\ 2 ^dxdt, 

where A is a Lagrangian multiplier enforcing |n| = 1, and the Oseen- Frank internal 
energy W is 

W(n, Vn) = T^MV • n) 2 + k 2 (n • V x n) 2 + A; 3 |n x (V x n)| 2 



^It is worth mentioning that the same kind of modification by a constant is required in the 
variational principle for the original HS equation in the periodic case, even without the /i-term. 
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where k± , k-z , and k% , are nonnegative physical parameters describing how the material 
responds to splay, twist, and bend, respectively. For the director field 

n(x, y, z, t) = cos ij>(t, x)e x + smip(t, x)e y , 

we get 

W(n, Vn) = i (ki sin 2 ip + fc 3 cos 2 2 . 

Writing 

c 2 (ip) = k\ sin 2 ip + k^ cos 2 0, 

it follows that 

W(n,Vn) = ±c 2 (^ 2 x . 
Hence, since | n £ | 2 = ip 2 , this leads to the action 

S = j j l -(tf-c 2 (^ 2 x )dxdt. 

(Note that the Lagrange multiplier A does not appear because n is already normalized 
to have length one.) The Euler-Lagrange equation, obtained in IIS, for the action S 
is 

iptt ~ c(ip)[c(ip)ip x } x = 0. 

At this point we deviate from the HS derivation and assume that there is a pre- 
ferred direction ipo of the director field. This can be related, e.g., to an exterior 
magnetic field, acting on dipoles and returning them to this direction. Such a return- 
ing force corresponds to the addition of the potential energy term — ^(ip — V>o) 2 to the 
Lagrangian density £ = \(4> 2 — c 2 (ip)ij) 2 ) . We also assume that there is retroactive 
action of other dipoles on a given director, i.e. that Hooke's elasticity coefficient 
k := fj,(ip — ipo) is proportional to the average deviation of all dipoles from ipo- other 
dipoles act on a given one, and the more other dipoles are rotated, the stronger the 
returning force is on the given dipole. 

Thus the new action is 

S = \f c 2 {^ 2 x - ^ - VoXV - V>o) 2 ) dxdt , 

where ipo is the direction about which the fluctuations occurU We need to choose V>o 
so that both Co := c(ipo), as well as c' := c'(ipo), are nonzero. 
The Euler-Lagrange equation now becomes 

(3.8) i> tt - c{$)W)$*\x + ~ ^)ty> ~ V>o) + \nM ~ V'o) 2 ) = 0. 

Now, repeating the derivation in [HS| with this modified equation, we introduce the 
variables 

9 = x — cot, r = et, 

and substitute the expansions 

Tp = Tpo + eipi(9,T) +e 2 ip 2 (0,r) + 0(e 3 ) and c{tp) = c + ec' Q ipi + 0(e 2 ) 
into (PT5|) . One finds 

(c + ec'o^i) (etpw + e 2 ip 2 e) 

. e 



(-c d e + ed T ) (etpi + e 2 ^ 2 ) - (co + ec' ipi) 



The added term does not preserve the periodicity of the angle if), but this is inessential since 
when deriving the asymptotic equation we only consider small fluctuations of ijj- 



8 BORIS KHESIN, JONATAN LENELLS, AND GERARD MISIOLEK 

+ (eVi + e 2 ^) J{eipi + e 2 ^ 2 )dx + (e^i + e 2 ^ 2 ) 2 dx + 0(e 3 ) = 0, 
where we used that 

d t = -c d e + ed T and d x = d e . 

Collecting terms of different order we see that both the O(l) and the 0(e) levels 
vanish, whereas the 0(e 2 ) terms give 

CoV>200 - ZcolpWr ~ C C [V>lV>16»] g ~ 0^289 ~ Co^lpllpiee + Ipl J ^\dx + \ J ^l dx = °- 

This is simplified to 

2 2c J 4c J 

Replacing r — ► t, 9 — ► x, and tpi — > u/c' , we arrive at 

"fas = -Jn 2 - uu 2a; + -—j-(un(u) + Jm(^ 2 ))- 
2 2c c 2 

This is the //HS equation in the form (|3.2p when 2e ^ e ; = 2. After differentiation this 
leads to the evolution equation (//HS) for ut xx . 



Then 



4. HAMILTONIAN FRAMEWORK AND CONSERVATION LAWS 
4.1. Bihamiltonian structure and the fiHS hierarchy. 

Theorem 4.1. The ijlHS equation is Hamiltonian with respect to two compatible 
Hamiltonian structures. 

Proof Let m := Au = fJ,(u) — u xx and define 

H\ = — J umdx = -(u,u), H2 = J (fi(u)u 2 + -uu 2 )dx. 

5H X 5H 2 , 2 , . , 1 2 

~5u =m ' -j^ = V{u ) + 2^(u)u- -u x -uu xx . 

Since A is self-adjoint, we have the relation = A^j^-. Hence 

(4.1) -El = u , = A' 1 (fi(u 2 ) + 2fj,(u)u - \u\ - uu xx ), 

dm dm 2 

so that the //HS equation 

^txx — 2[i{u?ju x -\- 2u x u xx -\- UU XXX i 

can be written both as 

R 6Hl 

dm 

and as 

R 6H2 

dm 

where the Hamiltonian operators are 

Bi = - (md x + 8 x m) , B2 = —d x A = d x . 

These operators are compatible as a straightforward but lengthy calculation shows. 
We will present their geometric meaning in the more general setting of the Virasoro 
algebra in Section [8l □ 
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Using this bihamiltonian structure it is well-known that one can (at least formally) 
obtain an infinite hierarchy of conservation laws, cf. [McK2j. If two functionals Fq 
and Fi satisfy 

dm dm 

we say that Fq raises to F\ or that F\ lowers to Fq; in symbols Fq /* F\ or F\ \ Fq. 
Starting from H\ and H2 the goal is to construct an infinite sequence of conservation 
laws such that 

• • • / #_i / Hq / Ht / H 2 / H 3 / ■ ■ ■ . 

Just like for the Camassa-Holm equation, we find that lowering of the H n 's is 
unobstructed for the /xHS equation, i.e. there are ff_i,iJ_2, . . . , such that \ 
H-2 \ • • • • Moreover, J3_i, i?-2, • • • are local functionals given by 

(4.2) H^^fm^dx, n = l,2,..., 

where &H ^~ 1 is obtained from S ^^ n via 

-= B i #2-^— = -^7= / -7= ( ) dx. 



bm 5m 2y/m J ^fm \ 5m 

Starting with 5H_\ = this provides a constructive approach for finding the 

£Z_ n 's, which is easily implemented in Mathematica. For example, we get 

5H-2 m xx 5m 2 
5m ~ 8m 5 /2 32m 7 / 2 ' 

5H- 3 1155m^. 23lmlm xx 2lm 2 xx 7m x m xxx m xxxx 



5m 1024m 13 / 2 128m n /2 64m 9 / 2 16m 9 / 2 16m 7 / 2 ' 
and the corresponding functionals H-2 and H-z are obtained from (14.2H . 

The situation for raising is more intricate due to the occurence of nonlocal func- 
tionals. Nevertheless, we can deduce the existence of expressions Q„[m] (possibly 
nonlocal in m) such that Qq = and 

Q n+ i = B^BiQn, n > 0. 

If there exist functionals H n [m] with = Q n , they constitute an infinite ladder of 
functionals with H\ /* H2 /* H2 f ••• ■ The existence of the H n is not evident in 
this approach due to the nonlocal nature of the Q n 's. However, it follows from the 
expansion of a Casimir for a family of compatible Poisson brackets in the approach 
of Section El cf. [KM] . 

In summary, the first few conservation laws in the hierarchy are 

H __!_/" ml ^ 5H^2 m ra 5m 2 



16 J m 5 / 2 ' 5m 8m 5 / 2 32m 7 / 2 ' 
8H-1 1 



H-i = J yfmdx, 



5m 2Jm 



Hq = / max, — — = 1, 

J 5m 

H\ = — f mudx, -— ^ = u, 

2 J 5m 

.2 , 1„A . SH 2_f.. «2\-l (..I. ,2\ , n..f..\.. 1 ..2 



H 2 = i n{u)u + -uu x J dx, ~§ m ~ = ( fI ~ d ^ {Ku) + 2fi(u)u--u x -uu a 
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m t 


=B 


m t 


=B 


m t 


=B 


m t 


=B 



XXX 



with corresponding commuting flows 

SH-2 = gfl-i = 
<5m 2 5m \^2i/m 

1— J — — U, 

om dm 

— = B 2 - — = -m x , 
dm om 

6Ht n 5H 2 
i- — = B 2 - — = -2mu x - m x u. 
dm om 

4.2. Relation with bi-variational principles. Recall that the ^HS equation can 
be obtained from two distinct variational principles, corresponding to the Lagrangian 
densities C and C introduced in Section El 

Theorem 4.2. The two variational formulations for the /J.HS equation formally cor- 
respond to the two Hamiltonian formulations of this equation with Hamiltonian func- 
tional Hi and H%. 

Proof The action is related to the Lagrangian by S = f Ldt and L = J Cdx. The 
first variational principle has the Lagrangian density of equation (|3.ip , 

C = u x ut + uu 2 + 2/j,(u)u 2 , 

where we skip the linear term. The momenta conjugate to the velocities ut are 

dC 

Hence the Hamiltonian density is 

7i = u x ut — C = —uu 2 x — 2fj,(u)u 2 , 

and so 

H = — J (uu 2 + 2fi(u)u 2 )dx, 
which indeed agress with the conservation law H 2 up to a factor of —2. 

In the second principle the Lagrangian density is 

- 1 1 

£ = 2 U l: + 2^(' u ) n + n ( Zt + UZx ) ■ 
The momenta conjugate to the velocities Ut,irt, and zt, respectively, are 

dC dC dC 

OUt OTXt OZt 

Consequently, the Hamiltonian density is 

H = z t ir - C = z t ir - Qu 2 . + ^n(u)u + ir(z t + uz x )\ . 
Using the equation of motion (|3.5p to substitute for irz x , this becomes 

H = ~-jU 2 x - uu xx + ^n(u)u. 
Therefore, the Hamiltonian is 

H = J Hdx = — J (u 2 + fi(u)u)dx, 
which is exactly H\. □ 
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Remark 4.3. (Lax Pair) The /jHS equation 

U txx - 2fl(u)u x + 2u x U xx + UUxxx = 0. 

can be viewed as the condition of compatibility between 

ipxx = Xmip, 

and 

where A is a spectral parameter and, as above, m = /u(it) — u xx . In other words, 
assuming that these equations are satisfied for all nonzero A G C, the relation iptxx = 
ip xxt holds if and only if u(t, x) satisfies the /iHS equation. 

5. The periodic Cauchy problem 

In this section we discuss well-posedness of the /jHS equation. We do not aim 
for the strongest possible results here. Instead, after establishing a basic local well- 
posedness theorem we proceed to show that while some solutions persist indefinitely 
other solutions develop singularities in finite time. 

Recall that the mean of any solution u = u(t, x) is conserved by the flow and hence 
the Cauchy problem for the /jHS equation can be recast in the form 

(5.1) u t + uu x + A- 1 d x (2fi(u)u + ^u x ^ =0 
and 

(5.2) u(0,x) = u (x), xeS 1 , iGl, 

where A = fj, — d\ is an isomorphism between i? s (S' 1 ) and H S ~ 2 (S 1 ) with inverse 
u = A~ l w given explicitly by equation (|2.3p . Conservation of the mean will play an 
important role later on. 

5.1. Local well-posedness and persistence result. In order to establish local 
well-posedness of the Cauchy problem (|5.ip - (|5.2p we adopt here the approach devel- 
oped in [M2] for the CH equation. Our main result is the following 



Theorem 5.1. If s > 3/2 then for any u a G i? s (5' 1 ) there exist T > and a unique 
solution u G C((-T,T),H S ) n C((-T,T), H 3 ' 1 ) of the initial value problem IQ) - 
h5.S\) which depends continuously on the initial data u Q . Furthermore, the solution 
persists as long as \\u{t,-)\\ c i stays bounded. 

First, however, we will study the associated Cauchy problem on the group V S (S 1 ) 
of Sobolev H s diffeomorphisms with s > 3/2. Given a solution u(t,x) of (|5.ip with 
initial data u Q we let t — > r](t, x) be the flow of u, that is 

(5.3) f]{t,x) = u(t,i](t,x)). 



Differentiating both sides of (|5.3p with respect to t we rewrite the original Cauchy 
problem (pUj)Hf5T2|) as 



(5.4) ^ = -|,4 1 d x (2fj,(rjori ^otj ! + !((?) 077 ^x) 2 ) j o 7] =:/(??, f}) 

(5.5) 77(0,0:) = a;, 77(0, x) = u (x) 
and prove 
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Theorem 5.2. Let s > 3/2. For any u Q G H S (S 1 ) there exist T > and a unique 
solution (77,77) £ C 1 ((— T, T), D S (S' 1 ) x .fP) 0/ i/ie initial value problem fa5.4\) - ll5.5\) 
which depends differentiably on the initial data u Q . 

An immediate consequence of Theorem 15.21 combined with the inverse function 
theorem for Banach spaces is the existence of the Riemannian exponential map in 
V s (S 1 ) whose geodesies were described in Theorem 1 2. II The corollary below permits 
a rigorous development of the associated infinite-dimensional Riemannian geometry 
and justifies the geometric constructions of Sections 2 and 7. 

Corollary 5.3. For any s > 3/2 the metric \2. 2\) has a well-defined exponential map 
exp id : U C Tid'D s (S 1 ) — > V s {S 1 ) which is a C 1 diffeomorphism from an open set U 
onto a neighbourhood of the identity in V s (S 1 ). 

Our first main task will be to show that the map (£, w) — ► (w, f(£,w)) defines a 
C 1 vector field on the Hilbert manifold V s (S 1 ) x i/T s (S' 1 ). The proof of this fact is 
similar to the corresponding proof for CH and we will refer the reader to [M2] for 
some of the details. 

Proof of Theorem 15.21 It is sufficient to show that / is continuously differentiable (in 
the Frechet sense) with respect to the variables £ and w in some neighbourhood of 
the point (id, 0) in T> s (S 1 ) x H S (S 1 ). The following technical result will be useful. 



Lemma 5.4. ([M2], Appendix 1) For s > 3/2 the composition map £ — ► w o £ with 
an H s function w and the inversion map £ — > ^ 1 are continuous from V s (S 1 ) to 
H S (S 1 ) and from V s (S 1 ) to itself respectively. Moreover 

(5.6) \\v°Z\\H'<C(l + \\t\\%.)\\w\\ B ' 
with C depending only on inf \d x £\ and sup\d x £\. 

In particular, using (|5.6j) we find that the right hand side of (|5,4j) maps (£,w) — > 
f(£,w) intotf^S 1 ). 

Next, rewrite / in the form 

(5.7) f^,w) = -A- 1 d x ^h^,w), 

where d x £ = dR^ o d x o dR^-i and Aj = dR^ o A^ 1 o dR^-i are the conjugations 
of d x and A^ 1 by a diffeomorphism £ (here denotes the right translation by £ in 
P S (S' 1 )) and where 

(5.8) h{t, w) = 2wJwo C X dx + \ 2 . 

Let s -> ( s be a smooth curve in V s (S 1 ) passing through the point £0 = £ i n the 
direction <9 s £ s | s= o = v. A formal computation gives 

dt(d X) t)(v) = d s (d x £ t )\ B=0 = [jjof^j.^^ 

and 

d^A' 1 )^) = dg(Aj*)\ a=0 = -Aj 1 d.(A i .)\ a ^ Q Aj 1 = -A-^voC'd^A]^- 1 , 

where [•, •] denotes the commutator of operators. Similarly from (|5.8p we find 
(5.9) 

d£hfc w -)(v) = 2{woC l J w ° ^d-xiv o -]^(d x {w o £~ 1 )) 2 d x {v o t' 1 )} o £ , 
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as well as 



dwh(£ w ){v) = 2< w o £ 1 / v o £ 1 dx + v o £ 1 / w o £ 1 < 
(5.10) ^ 7 

+ ia E (™or 1 R.( V or 1 )}°£. 



With the help of the above formulas we can compute the directional derivatives of / 
in (15.71) . namely 



(5.11) 

= {A- 1 d x (vo€- 1 d x (ho€- 1 ))-vo£- 1 d x A- 1 d x (ho£- 1 )}o£ 

-2^-% (w o r 1 / w o r 1 ^ o - l^w ° r 1 )) 2 ^ ° r 1 )) } ° e 



and 
(5.12) 



-2^A~ l d x (w o f" 1 J vo ^- x dx + v o f- 1 y io o f" 1 



+ ^(t»of 1 )^of 1 ))}oe 



Our strategy will be to show that both directional derivatives d^fr^ tW \ and d w fr^ w \ 
are bounded operators in L(H S ,H S ) that depend continuously on (£,w) in some 
neighbourhood of (id, 0). Using the identities 



(5.13) A 1 d x w(x) = (x — J w(x)dx — J w(y)dy + J J w(y)dydx 



and 



(5.14) A 1 d x w(x) = — w(x) + I w(x)dx 



which are easily obtained from (|2.3p and the fact that A 1 and d x commute, we 
find explicit expressions for the directional derivatives in (|5.1ip and (|5.12p . After a 
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lengthy, but straightforward, calculation we obtain 
(5.15) 



j)(v) =^4^x — Jwo£ 1 dx j vo C 1 9 x (w o ^ l )dx 

+ l( x ~\) I (^(« , °r 1 )) 2 ^(«°r 1 )^ 

/px p px 

w o £~ x d x (v o dx J wo^~ 1 dy — 2 J wo^~ 1 dx J v o ^~ 1 d x (w o dy 

-\f\d x (woC l )) 2 d x (voC l )dy 

+2 J wo^dx j j v o i~ 1 d x {w o dydx 

—2 / w o £ t ~ 1 d x (v o dx I I wo£~ 1 dydx 



o 



+\ j voC 1 {d x {woC 1 )fdx+ l - / /"^(wof'jj^of 1 )^ 



/ . , , .- - I . \ 2 , > , I 

X 







+2w o £ ^ o £ 1 y u; o £ 1 <ix — 2v o £ j w ° £ 
-\ V °Z~ X J {d x {woC l )fdx]oi 



2 



and similarly 



(5.16) 

+ ^) y ^x(w o C _1 )^a;( u ° £ _1 )^x — 2 y w o dx j v o dy 

/px px 
v o dx I w o ^ 1 dy — / d x (w o (,~ 1 )d x (v o £~ 1 )dy 
Jo Jo 

/p px p p px 

w o £ _1 dx / v o £ -1 dydx + 2 v o dx / / wo dydx 



By LemmaEH in order to show that v — > d^fr^ w \(v) is a bounded operator on 
it will be sufficient to estimate the sum 



Using Cauchy-Schwarz and the Sobolev lemma we can bound the first of the two 
terms above by 

< {\\w o c^Uh ° c^im + \\w o r 1 !!^) (Ik ° riifi + n« ° ricO 
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and, after differentiating (|5.15p in x and using the fact that is an algebra, 

estimate the second term by 

< \\d x (w o C l )\\ 2 L 2 \\d x {v o C 1 )^ + \\ d x(w o ^ 1 )||^ s -i \\d x (v o ^ _1 )|| Jf8 _ 1 

+ \\w ° C _1 llcx)||w ° C _1 ll-H" s ll u ° C 1 ^" + \\w o c^ll^Hd^ o ^ _1 )|| 2 j S -i 

Combining these and using the estimate (|5,6p of Lemma 15.41 we obtain 

II II ~ Clklllflbllir" > 

where the constant C depends only on the H s norms of d x £, and <9z£ _1 . An analogous 
argument yields a similar bound on d w f^^ in (|5.16p 

,mj)( W )||^ s CH^Hi/ 3 1 1 't-' 1 1 if s • 

We have thus shown that / is Gateaux differentiable near (id, 0). 

We next turn to continuity of the directional derivatives with respect to the vari- 
ables £ and w. Continuity in w follows from the fact that the dependence of both 
partials on this variable is polynomial. Therefore, it only remains to show that the 
norm of the difference \\d^f^ tW )(v) — df/fid,™)^)!!-^ can be made arbitrarily small 
whenever £ is close to id in T> s (S 1 ), uniformly in v, w G H s {S l ). Thus, our task is to 
estimate the sum 

\\ a ef(£,w)(v) ~ ^/(id llu )(u)|| L2 + \\d x (d{fe )VI )(v) - ^/(id.w)^))!!^-! • 

Using the explicit formulas in (15.150 together with the estimate of Lemma 15. 4} the 
algebra property of H^ 1 and the fact that for any r > we have a Sobolev embedding 
into Holder spaces ^(S 1 ) '—* C r ~ 1 > 2 (S 1 ) with a pointwise estimate 

\9(x)-g{y)\<\\9\\ H r\x-y\ r - 1 ' 2 (x,y G S 1 ) 

both terms in the expression above can be bounded by 

< C\\w\\%. \\v\\ H . (||r 1 - id\\ H . + Hr 1 - id||^ 3/2 ) , 

where C depends only on the H s norms of £ and its inverse. The arguments needed 
to prove continuity of the map £ — ► d w fr^ w \ are very similar. 

From these estimates we conclude that / defined by the right side of (|5.4|) is 
continuously differentiable in a neighbourhood of the point (id, 0). Applying the 
fundamental ODE theorem for Banach spaces we obtain Theorem 15.21 □ 



We are now ready to prove local well-posedness and persistence results for the /xHS 
equation. 

Proof of Theorem 15.11 Observe that directly from the flow equation f)5.3j) we have 
u = f] or] . Thererfore, since V s (S 1 ) is a topological group by Lemma 15.41 the first 
statement is a direct corollary of Theorem 15.21 

To prove the persistence result we employ Friedrichs' mollifiers J £ G OPS~°°, 
where < e < 1, and using (|5.ip we first derive an estimate for 




J £ u\\hs =-2(A s J e (ud x u),A s J E u) L 2 -4(A s A- 1 d x J £ ( f i(u)u),A s J £ u) 



- (A s A~ 1 d x J £ (d x uf,A s J e u) 
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where A s = (l — d 2 y^ 2 . The first term on the right hand side can be bounded as in 



\H S 



[T2] . see p. 115-116 

(A s J £ (ud x u) , A s J £ u) L 2 < \\u\\ c i\\u\\ 2 Hs . 
For the middle term we have 

(A 8 A~ 1 d x J e [fj l (u)u) , A s J £ u) L2 < \fi(u)\ WA'A-^JeuW^ ||A s J e u|| L2 < H^HI^. 
To estimate the last term we again use Cauchy-Schwarz 

{k s A~ l d x J e (d x u) 2 ,k s J e u) L2 < IIA- 1 ^^) 2 !^, \\u\\ H . 
and with the help of the identities (|5.13p and (|5.14[) we obtain 

WA-^dsufWx. ^ WA-^dsufWv + \\d x A- 1 d x (d x u) 2 \\ H .- 1 

< \\d x u\\h + II " (9 x u) 2 + J {d x ufdx\ 

< \\d x u\\ 2 L 2 + H^xuIIooIIuIIh* , 

where in the last step we used a standard Moser inequality for products of iP -1 
functions. 

Collecting these estimates and passing to the limit with e — > we arrive at the 
differential inequality 

d II m2 s 2 

j t \\u\\h° ~ IfIIc i IfIIh« . 

which combined with Gronwall's inequality implies that the H s norm of the solution 
u cannot blow up unless •) He 1 does. □ 

5.2. Global solutions with non-negative angular momentum density. Our 

next result shows the effect of the term 2fi(u)u x on the lifespan of solutions to the 
Cauchy problem. 

Theorem 5.5. Let s > 5/2. Assume that the function u a G H S (S 1 ) has non-zero 
mean and satisfies the condition 

(5.17) Au (x) := (n - d 2 x )u {x) > 0, x £ S 1 . 

Then the Cauchy problem 115. 1]) ~ {5J^) has a unique global solution u in the space 
C(R,H*)DC7 1 (R,H S - 1 ). 

Proof In Theorem 15. II for any such u Q we constructed a local solution u with desired 
regularity and defined up to some time T > 0. Furthermore, we showed that in order 
to extend it beyond T it suffices to find a constant K > such that \\u(t, Ollc 1 — K 
for all \t\ < T. 

Set m(t, x) = Au(t, x) and observe that from the identity (I5.13P we have 



m(t,y)dy 



d x u(t, x) = (x — — ^ J m(t,x)dx + J J m(t,y) dydx — J 

for any x E [0, 1], which consequently yields a bound 
(5-18) IIM*)lloo<||m(t)|| £ i. 

On the other hand, it is not difficult to check that at each point of the circle the 
solution u satisfies a local conservation law 

Au{t,n(t,x))(d x n{t,x)) 2 = Au {x) , 
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where rj is the geodesic in T> s {S 1 ) corresponding to u. In particular, if (|5.17p holds 
then m(t,x) > as long as u(t,x) exists, which in turn yields 

\\ mmL ,=ju { t, X ) dx = ju o{x)dx 

and the theorem follows. □ 

An analogous result holds if the assumption (|5.17p is replaced with Au {x) < 0. 
Notice that for the periodic Camassa-Holm equation the condition that does 
not properly change sign is necessary and sufficient for global existence, see |McK] . 

5.3. Breakdown of classical solutions. It is not difficult to find solutions which 
develop singularities in finite time. In fact, the conservation of the mean 

(5.19) J u(t,x)dx = J u (x) dx 

implies that sufficiently smooth solutions of /iHS corresponding to zero-mean initial 
data must satisfy the periodic HS equation 

(HS) Utxx ~\~ ^Uxllxx + UUxxx = i 

whose solutions inevitably break down. We recall here a simple proof of this facd 
Consider any smooth (non-constant) periodic solution u of the HS equation with 
smooth initial data u Q (see e.g. |Tij). Integrating (|HSp with respect to the space 
variable we obtain 

1 2 

Utx + + UU xx = a 

where 

a = — / u~ dx 
2 J 

is independent of time as is easily checked using periodicity and integration by parts. 
(Alternatively, one notices that a is a multiple of the i7 1 -energy of the solution, which 
is conserved for HS.) Since u' Q is continuous we can pick a point Xo where u' (x ) < 0. 
Then, along the characteristic t — > x(t) defined by x(0) = Xo we get an ordinary 
differential equation of the form 

— + -w 2 = a, where w(t) = u x (t, x(t)) 
whose explicit solution 

1 \ _ u 'o( x o)y/~ 2o — 2a tan (—t\J—2a/2) 
^ y^2a + u' (x ) tan (*V=2a/2) 

becomes unbounded as t f T cr u = /^ 2a arctan (\/—2a/\u' ( Xo )\) . 

Note, however, that blow-up of /iHS solutions can occur even if the mean of the 
initial data is not zero. More precisely, if u is a smooth solution of (|5.ip - (|5.2p then 
proceeding as above we arrive at the equation 

1 2 

u tx - 2fi(u)u + -u x + uu xx = r , 
where the constant of integration 

r = -2n(u) 2 - - / u 2 x d X 
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This argument was worked out by K. Gleason in her Senior Thesis [Gl] ; see also [Y] , 
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is again independent of time. Using Cauchy-Schwarz and conservation of the energy 
2H\{u) = /i(n) 2 + J u x dx and the mean (|5.19p . we then estimate 

u tx + \u 2 x + uu xx = r + 2fi(u)u 



2 



= 2fi(u )(^u(t,x) — J u{t,x)dx S j — — J {u' ) 2 dx 

< (2\n(u )\ - -IKH^^iki 



|L 2 

since the function u — fi(u) must vanish somewhere on the circle. The conservation 
law Hi ensures that solutions remain bounded as long as they are defined. Therefore 
the characteristics are defined as long as the solution exists. If we arrange for the 
right-hand side above to be non-positive then along any characteristic the function 
w(t) = u x (t,x(t)) will satisfy 

dw 1 9 

(5.20) - < --w* 

and so picking x Q as before with u' (x ) < we get 

w(t) < w(0) = u' (x ) < 

as long as it exists. Using this and solving (|5.20p we therefore obtain for any t > 
the estimate 

w(t)< 2<{Xo) 



2 + tu' {x ) ' 

from which conclude that w(t) becomes infinite at T cr u = 2/\u' (x )\. We summarize 
these results in 

Proposition 5.6. For any smooth initial data satisfying 4|^(-u )| < ||Uo||l2 there 
exists T cr u > such that the corresponding solution u(t, x) of the fiHS equation 
remains bounded for t < T cr a and satisfies ||ita;(i)||oo f 00 whenever t f T cr a. 

Thus, the /xHS equation is an example of a nonlinear PDE admitting solutions with 
finite as well as infinite lifespan^. As a simple illustration, the solution corresponding 
to the initial data u Q (x) = a+cos 2irx blows up whenever |a| < tt\^2/4 while persisting 
indefinitely if \a\ > 4ir 2 . 

It seems likely that Theorem 15.51 and Proposition 15.61 can be significantly strength- 
ened to give the full picture of well-posedness including a complete characterization 
of the blow-up mechanism (e.g. along the lines of |McK| or by studying the evolution 
of the extrema of u x as in [CE2j ). 

6. Traveling-wave solutions 

In order to find traveling-wave solutions of the //HS equation we first regard \i as 
a parameter independent of u. We look for solutions u(t, x) = ip(x — ct), traveling 
with speed c G R, of the equation 

(6.1) - u txx = -2/j,u x + 2u x u xx + uu xxx , 

where \x G R is a parameter. A first analysis will give us a list of all traveling waves of 
this equation. Subsequently, by keeping only the solutions (p for which the parameter 
[i equals the mean fJ,((p) = J tpdx of ip, the list will be narrowed down to traveling 
waves of the true /iHS equation. Since n(<p) is a conserved quantity, this procedure 
yields all traveling waves of ^HS. 



^Results of this type for the CH equation can be found e.g. in [CEj . 
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Substituting u(t, x) = ip(x — ct) into equation (|6.ip and integrating, we get 

1 2 a 

Ctfxx = -2fi(p + tptp xx + -ip x + -, 

for some constant d£l, Rewriting this as 

(6.2) & = -±W + {(<P-c) 2 ) X9 + a. 

we have an equation that makes sense for tp G (M) and it is natural to make the 
following definition: A function tp G fl^ c (R) is a traveling wave of equation HfiHS\) 
if it satisfies ( 16. 2p in distributional sense for some a G R. Applying the following 
lemma with v = tp — c we immediately deduce that any traveling-wave solution tp 
is smooth except possibly on the boundary of the set {x G R | tp{x) = c}. A 
general traveling- wave solution tp{x) is therefore smooth on the union of the (possibly 
countably infinite) number of open intervals that make up the complement of this 
boundary. 

Lemma 6.1. ([L]) Let pit) be a polynomial with real coefficients. Assume that v G 
H"* oc (R) satisfies 

(6.3) (v 2 )xx = v 2 x +piv) in P'(R). 

Then V s G C 3 (R) for k > 2 J . In particular, v is smooth away from the boundary of 
the set {x G R | v(x) = 0}. 

Within each interval where tp is smooth we may integrate (|6.2p once more to obtain 
the ODE 

2 -2^ 2 + a^ + o 

(6-4 = , 

c — p 

where b G R is another integration constant which may take on different values for 
different intervals. If the roots of the polynomial —2fip 2 + atp + b have nonzero 
imaginary part, there are no bounded solutions of (|6.4p . Thus, we may introduce real 
parameter m and M such that m < M and 

(6.5) -2/ip 2 + ap + b = 2(i(M -p)ip-m). 



The next theorem follows from an analysis of (|6.4p along the lines of [Lj. Note that 
(|6.2p is left unchanged under the substitutions 



(6.6) p. i— > — c/? i— > — (/?, c i— > — c, 

so it is enough to consider the case \x > 0. 

Theorem 6.2. Let p, > 0. For a given c£l, i/ie bounded traveling -wave solutions 
of speed c can be classified by the values of the parameters m, M G R as follows. 

(a) (Smooth periodic) If m < M < c, there is a smooth periodic traveling wave 
tpix — ct) of \6.1\) with min xg ]R p(x) = m and max^gjR p(x) = M. 

(b) (Cusped periodic) If m < c < M , there is a cusped periodic traveling wave 
tpix — ct) of W.l\) with minjjgR p{x) = m and max^gt p(x) = c. 

(c) (Anticusped periodic) If c < m < M , there is an anticusped periodic traveling 
wave pix — ct) of 116.1]) with min^gK p(x) = c and max ie i p(x) = m. 

(d) (Anticusped on the line) If c < m = M, there is an anticusped traveling wave 
pix - ct) of [61]) with min x .g R <p(x) = c and lim^^-too tp(x) = max^gK tp(x) = 
m = M. 
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(e) (Composite waves) Two wave-segments in (a) — (d) corresponding to param- 
eters (mi, Mi) respectively (1712, M2) are allowed to be glued together at their 
crests whenever mi + Mi = mi + Mi. Any composite wave formed by joining 
at most countably many wave-segments in this manner is a traveling wave of 

(/) (Plateaus) In the special case when mi + Mi = mi + Mi = c, the composite 
wave may also contain plateaus, i.e. intervals on which ip = c. 

6.1. Solitary waves. From Theorem 16.21 we infer the following list of solitary waves 
(i.e. traveling waves with decay to zero at infinity) of equation (|6.ip , 

• If fx > 0, then, for each c < 0, there is an anticusped solitary wave <p{x — ct) 
with min xe K <p(x) = c and lim^-toc, p{x) = 0. 

• If fx < 0, then, for each c > 0, there is a cusped solitary wave <p{x — ct) with 
lim^-too ip{x) = and max^R <p{x) = c. 

To verify that these are indeed anticusped respectively cusped solitary waves, we 
may determine their profiles directly. For example, choosing the case /i < so that 
c > 0, we employ (16.4p with integration constants a and b set to zero to get 



2 dip = x — xq. 



-2fup 

The integral can be computed to yield the formula 



(y/c — tp — ^/carctanh^y^— — —^j \ 



x — Xq. 



1 r 



U.5 



-1^ -1 -0.5 OJ 1 1.5 

Figure 1. Although this cusped solitary wave is a solution moving 
with speed c = 1 of equation (|6,ip when fj, = — 1, it is not a solution 
of the //HS equation due to the negative sign of //. 



Note that x — ► —00 as tp \ and that x — xq — > as cp f c. Near the crest point 
we get 

V2 



3cy/ — fx 



{c-<pf/ 2 + 0{(c-<pf' 2 ) = \x 



and so, near x = xq, 

c — ip(x) 



3c^—fi 

~7T 



2/3 



|x-x | 2/3 + O(|x-x | 4/3 ), 



showing that there is indeed a cusp at x = xq. 

Do any of these solitary waves satisfy the condition that \i equals J R ipdx, so that 
they are also solutions of the /iHS equation? The answer is no. Indeed, if /i > 0, 
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then since max^gR ip{x) = 0, the solitary wave has tp < everywhere, and so n((p) = 
f R ipdx < never equals \i. The case \i < is analogous. 

Theorem 6.3. The fiHS equation does not have solitary wave solutions on the line. 

There are also composite solitary wave solutions obtained by inserting cusped or 
anticusped periodic wave segments at the crests of these solitary waves according to 
(e) in Theorem 16.21 One might hope that a clever construction of a composite wave 
would yields a different outcome. However, since m + M = for the solitary wave 
and m < M by definition, one can see that any cusped wave segments joined to 
the solitary wave must have the same sign, and thus the sign of the integral J" R tpdx 
remains unchanged. 

This analysis shows that the fiHS equations admits no solitary waves, and hence 
no solitons, on the line. In the periodic setting for /iHS, however, we will see that 
traveling waves which are natural candidates for solitons do exist. 



6.2. Traveling waves on the circle. In view of the symmetry (16.6[) we restrict 
ourselves to the case fi > in this section. 
Consider equation ([67 



(6-7) <p x = 



C — if 

written in terms of m and M according to (|6.5p . By (a), (b), and (c) of Theorem 
there is a unique periodic traveling wave <p m ,M,u °f equation (|6.ip associated to each 
point (to, M, fj) such that (a) m < M < c, (b) m < c < M, or (c) c < to < M. By 
enforcing the constraint 

(6.8) mean,((p m) M,n) = A* 

we obtain the subset of these that are also solutions of the /uHS equation. From the 
way fi appears in (16.71) we infer that (see (16.121) below) 

(6.9) mean((p m , M ,fi) = -^=mean(<p m M 

The constraint (|6.8p therefore translates to mean(<p m ,M,i) = M 3 ^ 2 and is seen to admit 
a unique solution [i > exactly when Tnean(ip m 

,m,i) > 0. Whereas this condition 
fails to hold in the case of solitary waves as noted above, it is satisfied for large 
classes of periodic waves. Rather than finding the most general choices of m, M, and 
c for which it is fulfilled, we simply note that if the minimum of (p m Mi 1S positive 
then its mean is certainly also positive. Applying Theorem 16.21 with fj, = 1 we can 
therefore deduce existence of the following classes of traveling waves parametrized by 
m, M, c £ R. 

(a) (Smooth periodic) Whenever < m < M < c, there is a smooth periodic 
traveling wave ip(x—ct) of fiHS with min^gK tp{x) = to and max^gR tp(x) = M. 

(b) (Cusped periodic) Whenever < to < c < M, there is a cusped periodic 
traveling wave ip{x — ct) of ^HS with min xg R (p(x) = m and max^gR ip(x) = c. 

(c) (Anticusped periodic) Whenever < c < m < M, there is an anticusped peri- 
odic traveling wave ip(x—ct) of /xHS with min^gR (p(x) = c and max^gR tp(x) = 

TO. 

In the rest of this section we will derive explicit expressions (involving elliptic 
integrals) for the period and the mean of these waves. As a consequence of the 
analysis below, we will see that there are plenty of traveling waves of period one, i.e. 
on the circle S 1 = R/Z. 
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Theorem 6.4. For any c / the fiHS equation has a one-parameter family of smooth 
periodic traveling waves of period one and a one-parameter family of cusped periodic 
traveling waves of period one moving with speed c. 



6.2.1. Smooth periodic waves. From (|6.7p it follows that 



(610) z -* = 7%Lm-*Z-n)* f - 

Substituting ip = m + (M — m) sin 2 9, we find 



x — Xq 



H J 9o V c-m y (j, \ c-m J 



e 

which expresses x — xq in terms of the incomplete elliptic integral of the second kind 



E(d>\r) := / \A - r sin 2 6d0, -- < <> < 



7T , 7T 
- < 6 < — . 

2^2 



Now for smooth periodic waves the period is given by twice the integral from the 
minimum (fo = m to the maximum ipi = M, or, in terms of 6, from 9q = to Q\ = l|. 
We deduce that 



\2(c — m) ( M — ra 
(6.11) period(cp m ,M,iJ,) = 2W £ 



/U \ c — m 



where i?(r) := £7(7r/2|r) is the complete elliptic integral of the second kind. 
Similarly, we compute 



'2(c-m) /"^ , . /] M-m . . 

' / (7Ti + (Af-m)sm 2 (9)Wl sin 

Jd n V c-m 



V ^ Ve 
Then integration over a full period yields 
(6.12) 



. 2 2(c-m) ( , . ,„/Af-m\ , , r , r ,/M-m\\ 

= 3 A/ ^ (^(2(m + M) - C )S (^^^J + (c- AT)* 



where 

f7r/2 



/•7T/Z 

A"(r) := / (I - r sin 2 9)~ 1/2 d6, 
Jo 



is the complete elliptic integral of the first kind. Observe that the mean indeed scales 
with fj, in accordance with (|6,9p . 

Solving the constraint (|6,8p for // we see that the mean of the traveling wave cp of 
the //HS equation associated to the point (m, M) is 

, x (2 /- -(, , . s (M-m\ , , /M-m\\\ 2/3 

**>) = (3^^) (0(m + M) " c)* (— ) + (e - M)K (— ) ) ) 

and substituting this value of \i into ()6.1ip . we find that the corresponding period is 
(6.13) 

2 • W^fc^E 

period((p) 



(y^i ( (2(m + M ) - c)E (ff^) + (c - M)K (ff^))) 



1/3- 
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Figure 2. For each c / there are one-parameter families of both 
smooth and cusped period-one traveling waves of the /iHS equation. 



6.2.2. Cusped periodic waves. For cusped periodic waves (m < c < M) the period 
is given by twice the integral (|6.10p from ipo = m to ipi = c, i.e. from 9q = to 

We get 



arcsm 



M-m 



period(ip. 



m,M,jj, ) 



'2(c 



m 



E ( arcsin 



m 



M 



m 



M 



m 



m 



Similarly, formula (|6.12p for mean((p m; M,^) and formula f)6. 13|) for the period also 
hold in the cusped case with the replacements 



and 



where 



E 



K 



M 



m 



M 



m 



m 



m 



E l arcsm 



F ( arcsin 



m 



M 



m 



M-m 



m 



M 



m 



M 



777 



777 



m 



FW) := / (1 



rsm 2 d)-^ 2 de, 



is the incomplete analog of K. An analysis of the expression for the period in (|6.13D . 
together with its cusped analog, shows that the level set of period-one solutions looks 
qualitatively as in Figure In particular, we see that for each c > there are one- 
parameter families of both smooth and cusped period-one traveling-wave solutions of 
the nHS equation. This establishes Theorem! 



7. Geometry of the group of circle diffeomorphisms 

The objective of this section is to compute the curvature of V s (S 1 ) endowed with 
the fiHS metric (|2.2p . We establish a local formula for a connection compatible with 
the metric and use it to obtain an exact expression for the curvature tensor. In 
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particular, we prove that sectional curvatures in both vertical and horizontal two- 
dimensional directions on V s (S 1 ) are positive. 

7.1. Christoffel symbols and the covariant derivative. Define a symmetric bi- 
linear map T id (; •) on T id V s {S x ) ~ H s {S l ) by 

(7.1) T id (u,v) = —A~ x d x (ji(u)v + n(v)u + -u x v x ), u,v e H s (S l ), 

and extend it to a bilinear map r„(-, •) : T v V s (S r ) x T^P^S 1 ) -> T v V s {S l ) for any 
r/ £ P^S 1 ) using right invariance. 

Next, define the operator V locally in a chart at r\ by the formula 

(7.2) (VxY)fa) = £>Yfa) • X(r?) - ryFfr), X(r?)) 

for any vector fields X and Y on P S (S' 1 ). It turns out that V is the smooth Levi- 
Civita connection of the right-invariant /xrlS metric whose Christoffel symbols are 
given by the map V. The proof of the following proposition is straightforward (see 
for example [EM]). 

Proposition 7.1. The bilinear map V defined by is the unique torsion-free 

covariant derivative on V s (S 1 ) compatible with the right-invariant metric \2.2i) . 

7.2. Sectional curvature. In order to compute the sectional curvature of T> S (S 1 ) 
equipped with the fiHS metric (-, ■) defined in (|2.2p recall first that its curvature 
tensor R can be determined from the formula 

(7.3) R(X, Y)Z = V x VyZ - Vy V X Z - V [x ,y]Z 

for any vector fields X, Y and Z on V s (S 1 ). By right invariance it will be sufficient 
to compute R at the identity. We will write T(u,v) to denote the Christoffel symbol 
T id (u,v). 

Theorem 7.2. For any u,v £ T id V s (S l ) ~ H s (S l ), it holds that 

(R(u,v)v,u) ={T(u,v),T(u,v)) - (T(u,u),T(v,v)) - 3fi(u x v) 2 

(7.4) =/ u(u) 2 (/V) + M)) + (/V) + 

+ ii(u)fj,((vu x - uv x )v x ) + fj,(v)fj,((uv x - vu x )u x ) 

- 2(j,(u)fi(v)((j,(uv) + n(u x v x )) 

J u x v x dx^j + — J u x dx J v x dx — 3[i(u x v) 2 . 

Proof For any u, v and w £ T^ d T>(S l ) ~ H S (S 1 ) the local formula corresponding to 
the expression in (|T.3[) reads 

R(u, v)w =D 1 r id (w, u)v - DiT id (w, v)u 

+ r v (r id (w,v),u) -r id (r v (w,u),v) 

where D\ denotes differentiation of T v with respect to r\: DiTi d (w, u)v = 4:\ e =o^id+ev(w, it). 
We compute 

DiT id (w, u)v = -T id (w x v, u) - T id (w, u x v) + T id (w, u) x v, 

and so 

(R(u, v)v, u) =(—T(v x v, u) — T(v, u x v) + T(v x u, v) + T(v, v x u),u) 

+ (T(v,u) x v - T(v,v) x u + T(T(v,v),u) - T(T(v,u),v),u) . 
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Observe that 

F(u,v) = - ({uv) x - ad* v (u) - ad* u (v)), 

where ad*(u) = A~ 1 (2v x Au + vAu x ), see (|2.4p . Furthermore, note that the last four 
terms in the expression for (R(u, v)v, u) can be rewritten as 

(r(v, u) x - v(v, v) x u + r(r(«, v),u)- r(r(«, u),v),u) 

= (T(u,v),T(u,v)) - (T(u,u),T(v,v)) 
-{u x v, T(v, u)) + (uu x ,T(v, v)). 

Therefore 

{R(u, v)v, u) =(-T(v x v, u) - T(v, u x v) + T(v x u, v) + T(v, v x u),u) 
+ (T(u,v),T(u,v)) - (r(u,u),T(v,v)) 
- {u x v,T(v,u)) + (uu x ,T(v,v)). 

Now, using the formula for the Christoffel symbols in (|7.ip a straightforward com- 
putation gives 

(-T(v x v,u) - T(v,u x v) + T(v x u,v) + T(v,v x u),u) - (u x v,T(v,u)) + (uu x ,T(y,v)) 



This establishes the first curvature formula of Theorem 17.21 

To obtain the second equality in (|7.4R we simply substitute the expression for the 
Christoffel symbol T into 

(r(u, v),F(u, v)) — (T(u, u),T(v,v)) — 3( / u x vdx 



and integrate by parts to get 

/\ i 
(n(u)v + n(v)u + -u x v x )d x A~ d x (/j,(u)v + fi(v)u + -u x v x )dx 

+ J(2fi(u)u + ^u 2 x )d x A' 1 d x (2fi(v)v + ^vl)dx -3^J u x vdx\ . 

Using the identity d x A~ l d x f = A~ 1 d 2 = —f + fx(f) we simplify this expression and 
obtain the right-hand side of (|7.4p . □ 

From Theorem 17. 21 we can immediately deduce an expression for the sectional cur- 
vature K(u, v) of the plane spanned by two vectors u, v £ Ti^D s (S 1 ) . For this purpose 
we may assume, after taking linear combinations, that u and v are orthonormal with 
respect to the /^HS metric (•, •) and that fj,(v) = 0, i.e. 

Ktf+f «•*-!, „W=0, / = and /„ A «fa = 0. 

Js 1 Js 1 Js 1 

With these assumptions, Theorem 17.21 yields 

(7.5) K(u,v) = fi(u) 2 (n(v 2 ) + l) + n(u)fj,[(vu x — uv x )v x ) + ^ J u 2 dx — 3^(u x v) 2 . 

Proposition 7.3. For any orthonormal vectors u, v € TidD s {S l ) with fi(u) = /j,(v) = 
the sectional curvature K(u, v) of the plane spanned by u and v satisfies 

(7.6) K(u, v) = - - 3fi(u x v) 2 > 1 ( 1 - 1 ) > 0. 
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Proof In this case (|7.5p reduces to 
(7.7) K(u,v) =|-3^M 2 . 

By the Cauchy- Schwartz inequality it holds that 



J"' 

Since //(i>) = we also have 



v 2 dx < — =■ / i^cte 



4^2 y "x— 4 ^2' 

so that (USD follows from ([777]). □ 

Remark 7.4. The origin of the term — 3[i(u x v) 2 in (|7.6p is as follows. For any 

u, v G TidD s {S l ) decomposed as it = -u + ^(n) and i> = i; + //(/>) with //(it) = fj,(v) = 0, 
Theorem 17.21 yields 

(7.8) {R(u, v)v, u) = //-terms + {R HS (it, v)v, u) — 3[a(u x v) 2 , 

where '//-terms' stands for those terms containing either //(u) or fi{v ) (or both) as a 
factor, and 

(R HS (u, v)v, u) = — - ( I u x v x dx J +7 / u^dx / v^dx 

is the curvature of the H 1 metric (|2.ip corresponding to the HS equation, see |L2j . 
Then the middle term in (|7.8h corresponds to the curvature of the base P S (S' 1 )/S' 1 
contributing the constant term \ in (|7.6p . while the last term can be seen as the 
correction term in the general O'Neill formula |ON] 

(R M (X,Y)Y,X) = (R B (X,Y)Y,X) - ^||[X,y]^|| 2 , 

for a Riemannian submersion M — * B, where X,Y are vectors tangent to the base 
B, and X, Y are their horizontal lifts to M. Indeed, for horizontal vectors u and v 
the //-terms in (|7.8|) vanish so that since /j,(u x v) 2 = //([//, u]) 2 /4 the identification is 
clear. 

To summarize, let us decompose the tangent space at the identity TuV^S 1 ) = 
H © V in such a way that H consists of all mean zero functions on S 1 and V — R 
are the constants (so that u = u + /z(u) with //(it) = 0). Proposition 17.31 shows that 
the sectional curvature for any plane contained in (i.e. parallel to) the subspace H is 
strictly positive. Now turn to the planes perpendicular to H, i.e. planes containing 
constant functions, which constitute V. 

Proposition 7.5. Let v £ TidT> s (S ) be orthonormal to the constant function 1. 
Then the sectional curvature of the plane spanned by 1 and v is 



K(l, v) = J v 2 dx > 0. 



Proof Apply formula (|7,5p with u = 1 and v satisfying //(it) = 1, u x = 0, //(i>) = 0, 
and f sl v 2 dx = 1. □ 
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Remark 7.6. Applying Proposition 17.51 to the functions {vj.(x) = j^ kx }kLi, we 
obtain a sequence of two-dimensional subspaces whose sectional curvatures K(l, Vk) = 
p- approach zero as k — ► oo. 

Remark 7.7. The group V s (S 1 ) can be thought of as a solid torus, with a cross- 
section isomorphic to T> S (S 1 )/S 1 . As we mentioned in the introduction, there is a 
natural Riemannian submersion of the metric (•,•) on T> s (S 1 ) to the H 1 -metric on 
V S (S 1 )/S 1 defined in (|2.ip . The latter is isometric to a (contractible) piece U C S°° 
of the unit sphere S°° in L 2 (S l ), see |L2] I^ This set U has a finite diameter, and any 
geodesic on it reaches the boundary of U before acquiring a conjugate point. 

This geometric picture sheds light on the (non)existence of global solutions proved 
in the section on well-posedness. Since the group T> S (S 1 ) with the metric (•,•) is a 
solid torus whose cross-section is a slight deformation of the open subset U (of finite 
diameter), it is natural to expect that geodesies starting almost parallel to the cross- 
section (i.e. with small values niuj] would run into the boundary of the solid torus 
(and only deviate slightly from the HS geodesies, which hit the boundary in short 
time). On the other hand, geodesies starting 'perpendicular' to the base would wind 
around the torus indefinitely and correspond to global solutions. This is exactly the 
picture which was made precise in Section [5j 



8. //HS AS A BI-HAMILTONIAN SYSTEM ON THE VlRASORO GROUP 

In Section [2] we showed that the //HS equation is the equation of the geodesic 
flow of the right-invariant metric (|2.2[) on the diffeomorphism group T> s (S 1 ), see 
Theorem 12.11 In other words, the /iHS equation is the Euler equation on the Lie 
algebra!] Oect = T^V^S 1 ). By using the inner product (12. ip to identify Oect with its 
dual algebra Dect*, the equation can be viewed as Hamiltonian with respect to the 
canonical Lie-Poisson structure and the Hamiltonian H±. 

In this section we shall carry out this construction on the Virasoro group, which 
is the universal central extension of P(S' 1 ) or, equivalently, on the corresponding Lie 
algebra called the Virasoro algebra. As with the KdV and CH equations the Virasoro 
group and algebra provide a natural setting for a geometric interpretation of the bi- 
Hamiltonian structure of the //HS equation discussed in Section [H We will see that 
passing to the Virasoro algebra introduces an extra term 2ku x . The corresponding 
Euler equation reads 

(8.1) - utxx = -2fi(u)u x + 2ku x + 2u x u xx + uu xxx (k G R). 

The situation is analogous to that of the CH equation which on the Virasoro algebra 
becomes (cf. [KM] ) 

Ut ~ Utxx + 2ku x + 3uu x = 2u x u xx + uu xxx (k G R). 

In this case the parameter k is related to the critical shallow water speed, see [J]E3 



^One can show, however, that there is no isometric embedding of the group V s (S 1 ) endowed with 
the /xHS metric to S°° x S 1 , equipped with the direct product metric. 

Strictly speaking, Ti < j© s (5' 1 ) may be considered a Lie algebra only when s = oo. The reader 
can either make this assumption throughout the section or else take the Sobolev index s sufficiently 
large for all constructions to be rigorously justified. 

^One should also mention that in the case of the CH equation there is an important difference 
in regularity properties of the associated Riemannian exponential maps on the Virasoro and the 
diffeomorphism groups, see [CKKT] . 
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Recall that the Virasoro algebra is the vector space trit = Dect ® R endowed with 
the following commutator between pairs 



where v, w £ dect and &,cGl. Define an inner product on this algebra by 
(8.2) {(v,b),(w,c)) = (v,w) +bc 

where (•, •) is the inner product in (|2.ip . and use it to identify the dual trit* with trit 
by the pairing 

{ (m,a),(v,b)y=Jm vdx + ab 

where (u, a) i— > (Au, a) = (m, a) is the corresponding inertia operator. 

Let h, f : trit* — * R be arbitrary smooth functions. The dual space trit* carries the 
canonical Lie-Poisson bracket 

m {*./}(™,.) = <(™,.),[(^^),(^,g)])-. 

On the other hand, given any point (mo, ao) in t>it*, we can also define 

(8.4) {h,fUm,a) = ((m ,a ), \(^),(fMY 

\ l\om oa/ \dm oa/\l 

called a frozen (or constant) Poisson bracket. It is not difficult to check that both 

brackets form a Poisson pair in the sense that a linear combination of these brackets 

is again a Poisson bracket. 

Our main result in this section is 



Theorem 8.1. For any k £ R the /J.HS equation ( Iff. 1\) is Hamiltonian with respect to 
the two Poisson structures §8.3\) and j8.4\) on trit* where the second bracket is frozen 
at the point (mo, ao) = (0, 1)- 

Proof From the relations 

( a d* Vtb )(m,a),(w,c))* = ((m,a),[(v,b),(w,c)})* 

= ((m,a),(v x w - vw x , J v x w xx dxj^ 

holding for any elements (m, a) G trit* and (v, b), (w, c) G Dtt, we deduce an explicit 
form of the coadjoint action on trit* 

(8.5) ad* v b j(m, a) = (m x v + 2mv x + av xxx , 0). 

Since the Euler equation on trit* is of the form (m, a)t = —Q-d^ u \(m, a) we obtain the 
system 

{nit = —m x u — 2mu x — au xxx 
(H = 0. 

The second of the above equations immediately implies that a must be constant and 
using m = Au the first equation can be rewritten as 

^txx — 2 fjj(uju x -\- t 2u x u xx -\- uu xxx au xxx . 

Substituting u — * u + a we recover the /iHS equation in the more familiar form 
(18. ID with the constant k = —a. Since the Euler equation on trit* is automatically 
Hamiltonian with respect to the Lie-Poisson structure (|8.3h and with 

'2 X 



H l{ m,a) = \^iuf + Juld X + 
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as its Hamiltonian function, the first part of the theorem follows. 

To complete the proof, observe that the corresponding Hamiltonian equation for 
the frozen Poisson bracket (j8.4j) and any function h : trit* — > M reads 



(m,a)t = -ad*,sh_ shJm ,a ), 

\ Sm ' 5a J 



which, with the help of the formula (|8.5|) for the coadjoint action, gives 

Sh n ( Sh \ / Sh \ 

m t = -m 0x - 2m — - a — 

dm \dm/x \ dm/ xxx 

(h = 0. 



(8.6) 

Next, choose the Hamiltonian function h to be 



ft(m,.)-/(*«y + i«4 + *n»)* 



and compute its variational derivative 

= A' 1 (^i{u 2 ) + 2fi(u)u - ~v%, - uu xx + km 

cf. equation (14. ip . Substituting this expression into (]8.6p and picking (mo,ao) = 
(0, —1) gives 

J mt = c^A -1 {^{u 2 ) + 2/j,(u)u — \u x — uu xx + fcm) , 
\ a t = 0. 

Since c^A -1 = — 9 X and m t = —ut X x, the first of the above equations becomes 

^txx — 2fi(u^u x -\- 2u x ii xx -\- wu xxx ~t~ ku xxx 
and a substitution u — > u — k gives equation (|8.1[) . □ 



Remark 8.2. According to the classification given in [KM], on the Virasoro dual 
there are generically (i.e. at each point of the 'highest-dimensional' Virasoro coadjoint 
orbits) only three types of Poisson pairs consisting of the Lie-Poisson and frozen 
brackets. Each such pair is represented by one of the three integrable equations: 
KdV, CH and HS. On the other hand, Theorem 18.11 shows that the /iHS equation has 
the same bihamiltonian structure as KdV and corresponds to the same Poisson pair. 
How does this new integrable equation fit in the classification of [KM]? 

Observe that a Poisson pair determines a bihamiltonian equation only up to a 
choice of a Casimir (see Corollary 6.4 of [KM]) and this choice is different for the 
^HS and the KdV equations. In fact, it is possible to choose Casimirs at each step 
of the Lenard-Magri scheme. In the case of //HS and KdV we start with the same 
Casimir J udx (cf. Hq in Section 14. ip but the subsequent choices of Casimirs are 
different, which results in two different integrable equations albeit having the same 
bihamiltonian structure. 
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